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Abstract. The ringdown phase of a black hole formed from the merger of two orbiting 
black holes is described by means of the close-limit (CL) approximation starting from 
second-post-Newtonian (2PN) initial conditions. The 2PN metric of point-particle 
binaries is formally expanded in CL form and identified with that of a perturbed 
Schwarzschild black hole. The multipolar coefficients describing the even-parity (or 
polar) and odd-parity (axial) components of the linear perturbation consistently satisfy 
the 2PN-accuratc pcrturbativc field equations. We use these coefficients to build initial 
conditions for the Regge- Wheeler and Zerilli wave equations, which we then evolve 
numerically. The ringdown waveform is obtained in two cases: head-on collision with 
zero-angular momentum, composed only of even modes, and circular orbits, for which 
both even and odd modes contribute. In a separate work, this formalism is applied to 
the study of the gravitational recoil produced during the ringdown phase of coalescing 
binary black holes. 
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1. Introduction 

Post-Newtonian (PN) methods liave proved to be extremely efficient in describing the 
inspiral phase of compact binary systems, up to about the location of the innermost 
circular orbit (ICO). The PN inspiral signal has been developed up to 3.5PN ordeijf] 
for the orbital phase evolution [H |2] and up to 3PN order P HJ |5] in the amphtude 
waveform (see [6] for a review). On the other hand, recent advances in numerical 
calculations of binary black holes [Tj El [9] have provided a very accurate description of 
the subsequent merger and ringdown phases, say, from the ICO on. The comparison 
of the numerical-relativity and PN results is a crucial task that has been successfully 
achieved [TOl HH [121 [131 HI]- Their matching is currently under way [I5] and should 
yield a complete and very accurate solution of the problem of binary coalescence. 

Nevertheless, analytic and/or semi-analytic methods are still very useful for gaining 
more physical understanding of the relaxation of binary black holes toward their final 
equilibrium state (see e.g. [12] for a recent example). Of particular importance is the 
close-limit (CL) approximation method, pioneered by Price & PuUin |T7j and Abrahams 
& Price [IH]. The CL approximation permits the description of the last stage of 
evolution of a black hole binary, when the two black holes are close enough that they are 
surrounded by a common horizon, as a perturbation of a single (Schwarzschild or Kerr) 
black hole. Recent revisits of the CL approach made use of numerically generated initial 
data \l9l , and Bowen- York-type initial conditions [20l [21] . Alternative analytic or semi- 
analytic schemes for dealing with the same problem are based on the effective-one-body 
approach [22] (see [23] for a recent review). 

In the present paper we shall implement the CL approximation starting from post- 
Newtonian initial conditions, appropriate for the initial inspiral phase of binary black 
holes. A physical motivation is that the results of numerical relativity [71 [H [9] show that 
the pulse of radiation coming from the merger phase is very short and seems to connect 
smoothly to the previous inspiral and subsequent ringdown phases. It is thus reasonable 
to expect that PN initial conditions starting the CL evolution should essentially yield 
the right physics for the ringdown phase. The application of this formalism to the 
computation of the gravitational recoil effect or "kick" occuring during the ringdown 
phase will be presented in a separate work [21] . 

Let us outline the method. We conveniently distinguish several dimensionless ratios 
to describe a compact binary system. First, we introduce the "post-Minkowskian" (PM) 
ratio, measuring the internal gravity responsible for the dynamics of the system, and 
defined by 

GM . , 

£PM ~ ^3 , (1-1) 

where ri2 is the typical distance between the two compact bodies, and M = mi + 1712 

I As usual the nPN order refers either to the terms ~ in the equations of motion, with respect 

to the usual Newtonian acceleration, or in the radiation field, relatively to the standard quadrupolar 
waveform. 
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is the sum of their masses. Second, a post-Newtonian expansion will essentially be an 
expansion in powers of the a priori distinct slowness parameter 

ePN~%, (1.2) 

where Vi2 is the typical value of the orbital relative velocity. Recall that the PN 
expansion is only valid in the near zone defined by r -C A, where A ~ ?"i2/a/^pm is 
the typical wavelength of the emitted gravitational waves, and r the distance from the 
field point to, say, the center of mass of the binary. 

For a binary system moving on a circular orbit the two parameters epN and epu 
are comparable, epN ~ Spu- In this case, if we limit the PN expansion to a few terms, 
we need epN ^ 1 hence ri2 ^ GM/c^. However, it is often better to view ep^ and epu 
as independent parameters because if the binary system is moving on a highly eccentric 
bound orbit with eccentricity e < 1, the PN parameter can be much smaller that the 
PM one, since at the apoapsis of the orbit we have £pn ~ (1 — e) epu <^ £pm- And, for 
an unbound orbit with eccentricity e ^ 1, we would have epN ~ (1 + e) ^pm ^ ^pm at 
the periapsis. 

On the other hand, the close-limit approximation consists of an expansion in powers 
of the dimensionless ratio considered smal^ 

BCL^—- (1-3) 

r 

This expansion can formally be viewed as an expansion when the size of the source tends 
to zero, or multipolar expansion. Therefore, if we limit the expansion to a few terms, we 
need Eql ^ 1 and the CL approximation is expected to be valid in the domain r ^ 
(like a multipole expansion). 

Clearly the PN and CL approximations that we intend to employ simultaneously 
have disconnected domains of validity. Indeed, the CL describes a slightly distorted 
black hole such that ri2 > GM/c^, so that for circular orbits epN ^ 1, which makes 
the near zone very small; in other words, the PN metric will only be valid very close to 
the source while the CL approximation requires r ^ ri2. Despite such apparent clash, 
we shall be inspired by the method of matched asymptotic expansions ^3], which in 
principle allows one to get an analytic expression valid in the entire domain ^ r < +oo. 
Of course, this method is based on the existence of an overlapping zone, where the two 
asymptotic expansions are simultaneously valid and can be matched together. But in 
the present context there is no such thing as a overlapping zone. Hence our use of the 
theory of matched asymptotic expansions to relate PN and CL approximations can at 
best be only formal. 

Starting from the PN metric, already in the form of an expansion in powers of £pN, 
we shall restrict ourselves to the terms linear in epu (i-e., essentially, linear in G). This 
is to be consistent later with the use of a linear black hole perturbation. Then, each of 
the coefficients of the PN metric will be expanded in powers of Ecl, which will enable 

§ In the works [17l[T8], the CL parameter is defined as c'^ri2/{GM). More recently, Sopuerta et al. 
[20] adopted the definition (1.3). In the formal limit ri2 — >■ 0, these two definitions are equivalent. 
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us to identify the Schwarzschild background metric (up to terms of order Epy^) and the 
perturbation h^^, of that background. Thus the perturbation will appear as a double 
expansion series of the type [cf. the explicit results (2.12a)-(2.12c) 

h,. = epMj211 ^ 41' + 0{el^) , (1.4) 

where n refers to the post-Newtonian order and k can be viewed as the multipolar 
order of the expansion|j]] In principle, one could perform the expansions in the opposite 
way, i.e. expanding first in powers of ecL, and then in powers of £pn. In the method 
of matched asymptotic expansions the result should be the same, i.e. term by term 
identical in the double expansion series. This would however require first the knowledge 
of the black hole perturbation metric in the CL approximation; such metric can only be 
computed numerically. 

In the present paper, we shall implement the expansion (1.4), limiting ourselves to 
second post-Newtonian order. The reason is that the metric is needed in closed analytic 
form for any field point in the near zone, and that the 3PN metric is currently not 
known for any field point; only the 3PN metric when regularized at the very location of 
the particles is known One of our aims is the study reported in the separate work 
[21] of the gravitational recoil effect. The recoil is the reaction of the binary system to 
the linear momentum carried away by the gravitational waves, and results at leading 
order from the interaction between the 1 = 2 and the £ = 3 modes, where i is the 
azimuthal number of the decomposition of the black hole perturbation onto tensorial 
spherical harmonics; we shall thus push the CL expansion up to at least octupolar order, 
i.e. /c ^ 3, to ensure that the modes 1 = 2 and 3 are both taken into account. 

The present approach will be limited to the case of a slowly spinning black hole. The 
initial orbital angular momentum of the binary system, which is constant and supposed 
to give the spin of the final black hole, is considered as part of the perturbation of 
a Schwarzschild black hole, and is necessarily small. However we have learned from 
numerical calculations that the final black hole produced by coalescence is a rapidly 
spinning Kerr black hole [Zl [HI E] . The ringdown waveform that we shall compute in this 
paper will be that of a perturbed Schwarzschild black hole, and hence the quasi-normal 
mode frequencies will not include the effect of the black hole spin. To remedy this 
problem and get better agreement with numerical relativity would necessitate similar 
calculations using a Kerr black hole background. 

The remainder of this paper is organized as follows: In Sec. [2] we consider the 
2PN metric of two compact bodies at first post-Minkowskian order (to be consistent 
with first-order black hole perturbations) and formally re-expand it in the CL form. In 
Sec. [3] we first give a short recap of the theory of linear perturbations of a Schwarzschild 
black hole, and then use this formalism to identify the perturbation associated with 
the previously CL-expanded 2PN metric. In Sec. [4| we verify that the field equations 

II Our convention is that fc represents the power of the separation ryi hi the CL expansion, taking into 
account the inverse power of ryi hidden in the PM indicator epM in front of (1.4). 
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for this metric are satisfied. Using the CL-expanded 2PN metric as initial data, we 
numerically evolve the Regge- Wheeler and Zerilli functions in Sec. [5| and present the 
resulting waveforms generated during the ringdown phase of coalescing black holes, for 
both even and odd-parity perturbations. Finally we conclude in Sec. |6| Some necessary 



details on black hole perturbation theory are relagated to Appendix A 



2. The 2PN metric in close-limit form 

2.1. The 2PN metric for two point masses 

In the present paper we shall solve numerically the Regge- Wheeler and Zerilli wave 



equations [see Eqs. (5.3) below] starting from post-Newtonian initial conditions. Thus 
we assume that the initial metric at the end of the inspiral phase is given by the standard 
PN metric generated by two point masses mi and m2 modelling two non-spinning black 
holes. We adopt the 2PN precision because the 3PN metric in the near zone is not 
known in the "bulk", i.e. for any field point outside the position of the particles. 

Our calculation starts with the post-Newtonian metric G™ written in a Cartesian 



harmonic coordinate system, and given as j27|^ 



PN _ 2Gmi 1 

*-^00 — ~ -L H 2 ^ ~^ 



ri 

2 



+ G^mims + ^ 



rir2 1r\2 Ir^r^i. 2r2ri2 



AG mim2 



+ ^ , 2 (^12^12) + 1^2 + 0{c-^) , (2.1a) 



Gnii ■ , 1 



'^Oi ~ ^ ^3^. ^1 + ^5 



G'^ml G'^mim2 I ^ , .^r \ 

^1 ( ;:2"^(^i^i) + — g-2 ( -16(ni2t;i) + 12(ni2f2) 

-16(n2fi) + 12(n2f2' 



/ ri 1 1 

+ n\^G'^mim2[ -Q{ni2Vi2)— - -iiniVi)— + l2{niVi) — 

-l6K..)i+4(n,.,)i(i + i- 



+1' 



Grrii , , r,x G^mi /Sri 2r2 

— ^ {2{niv,f - Avl) + — ^ + G'm,m2 [-^ - 

'1 '1 V'l2 '12 



^ Greek indiees take space-time values 0,1,2,3. The indices fijV,... indicate Cartesian coordinates 
a;'* — {ct, X, y, z}, while a, l3, . . . refer to spherical coordinates — {ct, r, 9, if}. Latin indices i, j, . ■ ■ 
take spatial values 1, 2, 3. Bold-face notation is often used to denote ordinary spatial vectors, x — {x*}. 
The two black holes are often labeled by A = 1,2. Parentheses around indices are used to indicate 
symmetrisation, e.g. U^^V^'' — \ (U^V^ + t/^V'). The usual (Euclidean) scalar product between two 



3-vectors U and V is denoted {UV), e.g. {nivi) = rii • vi. To the terms given below in Eqs. (2.1a 



(2.1c I, we have to add those ones corresponding to the relabeling 1 o 2 (with the obvious exception of 



the Minkowski metric which should not be counted twice). 
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-G'^mim2 



+ 



riri2 r2ri2 risS" 



5'^ + 



+1^2 
2Gmi 



0{c 
1 



(2.16) 



+4 



Gmi 
ri 



v\v{ + 



V ^1 



2 



rir2 



12 



n\n\ — 4(7 mim2n\2n 



12 



1 

52 



2riri2 
1 



+ 



^125' 



^25 



+ 



4:G'^mim2 



g2 



n]^ 712 + ^^1 '^12 



G'^mim2 
&rl2 



-l^{riY2VY2)S'^ - 6(ni2t;i2)^i2 



12 



f2.1c) 



+1^2 + Oic-") . 

The trajectory of the A^^ black hole is denoted yA and its ordinary velocity is 
va = dyyi/dt, where t = x^/c is the harmonic-coordinate time. The relative velocity 
is denoted V12 = vi — V2. The Euclidean distance between the black hole A and any 
field point is = |x — y^l- The unit vector pointing from A to the field point is 
= (x — jAj/rA, and the unit direction from body 2 to body 1 reads ni2 = yi2/''"i2, 
where yi2 = yi — y2 and the binary's separation is denoted ri2 = |yi — y2|- Several 
terms in (2.1a)-(2.1c) involve the particular combination S = ri -\- r2 + r\2- AH these 



conventions can be visualized in Fig. [T] 

In the following we shall restrict ourselves to those terms in the full 2PN metric 



(2.1a)-(2.1c) which are linear in the parameter epM given by (1.1), or equivalently in 



the gravitational constant G^\ Indeed, our work will be based on the theory of first- 
order perturbations of a Schwarzschild black hole, and the corresponding terms in the 
PN framework will necessarily have to involve only linear powers of G to be consistent. 
Higher powers of G in the PN metric will correspond to higher-order perturbation theory. 
Since the part of the 2PN metric which is linear in G is obviously a solution of the 
Einstein field equations at 2PN order [up to terms (^(G^)], the multipolar coefficients 
describing the perturbation in the CL approximation will satisfy the perturbative 
Einstein field equations, as checked in Sec. |4j This restriction to the terms linear 
in G in the PN metric appears therefore as necessary; the price we pay is that our 
initial conditions will not contain the full information encoded into the 2PN metric: for 



instance, all the terms involving S in (2.1a)-(2.1c) disappear. To include meaningfully 



the complete 2PN metric would require using the theory of second-order perturbations 
of a Schwarzschild black hole |28j . 

Neglecting the non-linear terms in G we end up with a comparatively much simpler 

+ From now on it will be simpler to forget about the dimensionless estimates epM , £PN and £cl defined 
in the Introduction for pedagogical reasons. We shall replace them by the dimensionful but more 
obvious constants G and c~^, and parameter r\i, respectively. 



Close-limit approximation with post- Newtonian initial conditions 



r 




X 



n 



m 



Figure 1. Two black holes of Schwarzschild masses mi^2 in a Cartesian coordinate 
system = {ct, x, y, z\ centered on the binary's center of mass, or the associated 
spherical coordinate system — {ci,r, 0, (^}. 



metric, reading 

2Gmi 



^PN 
'-^00 



/-yPN 



-1 



1 + 



AGmi 



2vf 
1 



2 



nivi 



1 ^2 + 0{G^,c 



1^2 + 0{G^ 



1 + 



2Gmi 



Gnii 2 



6'^ + ^^v\v{ + 1^2 + 0{G\ c~^) 



(2.2a) 
(2.26) 
(2.2c) 



At this stage we could proceed with the CL expansion to identify the Schwarzschild 
background metric and the perturbation. However the PN metric is in harmonic 
coordinates so we would obtain the Schwarzschild metric in harmonic coordinates; 
this is not convenient because the perturbation formalism is usually defined in 
standard Schwarzschild-Droste coordinates. We shall thus perform a suitable coordinate 
transformation such that after expanding the metric in the CL approximation we obtain 
directly the Schwarzschild background metric in Schwarzschild coordinates. 

Since we are working at linear order in G it is sufficient to perform a linear gauge 
transformation at order G, say 5x^ = where the gauge vector is ^'^ = 0{G). Note 



* The remainder 0{G^,c~^) includes all terms which are at least of order or of order or both. 
Thus for instance the "radiation-reaction" terms present at order c 
this remainder because they are also of order G^. 



m Ctqq 



and GfJ^ are included in 
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that this gauge transformation is defined with respect to the Minkowski background. 
Later, when using black-hole perturbations, we shall perform a gauge transformation 
with respect to the Schwarzschild background. A suitable gauge transformation which 
fulfills our purpose is then 



e 



Gmi 



Gm2 



(^2^2) 



Gmi ■ Gm2 ,■ 



(2.3a) 
(2.36) 



Such gauge transformation does not satisfy the harmonic gauge condition = 0, 
therefore our new coordinates will not be harmonic. Under this gauge transformation 
we have 



G'^-[x] 



where 



f]^u^^, and the new metric components read now 



PN 



9^ 



-1 + 



2Gmi 



1 + 



3^? - 2(^1^0 



+ 1 ^2 + 0{G\c- 



AGmi 



c^ri 



[nivi, 



v{ + 1^2 + OiG^,c- 



6'^ 



2Gmi 



n\n\ + 



2v\v{ - -{n^v^fd^^ 



1 ^ 2 + C(G 



^c-«) 



(2.4) 

(2.5a) 
(2.56) 
(2.5c) 



We shall start our perturbative CL setup from that PN metric. 

It is important to comment on the post-Newtonian counting we are adopting for 



the 2PN metric in (2.5a)-(2.5c). The accuracy of this metric is really 2PN only for 



the geodesic motion of "photons" rather than of massive particles. Indeed, we did not 
include here the term of order 0{c~^) in the 00 component of the metric, although 
it is known from [27]. This term would be needed for describing the 2PN motion of 
massive bodies. It will turn out to be essential to expand both g^^ and f?™ at the 
same post-Newtonian order — namely up to ~ for both components in (2.5a) and 



(2.5c) — because only then can we be consistent with the linear black hole perturbation. 



Physically this results from the fact that the CL approximation is assuming that the 



metric is a small deformation of that of a black hole, therefore when expanding (2.5a)- 
(2.5c) in the CL form we shall have ri2 > GM/c^, so the orbital velocities are going to 



be very relativistic, i.e. \wa\ ^ c, and thus the 5^00 , 9ij and g™ components should give 
similar contributions to the line element ds^, and are therefore to be expanded up to 
the same PN order. This is thus similar to the accuracy needed for the geodesic motion 
of a photon where all the metric coefficients (7™, (7™ and g™ are to be given with the 
same PN accuracy. On the other hand, we shall see later that it is very important to 
include the gravitomagnetic field g™ up to order c~^, because it will play a crucial role 
in generating the odd-parity or axial contributions to the perturbations. 
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2.2. Close-limit expansion of the 2PN metric 



In the CL approximation, we expand the 2PN metric (2.5a)-(2.5c) in powers of the 
parameter (1.3), or equivalently of the relative distance ri2 = |yi — y2| when ri2 — t- 0. 



To do so we shall first express all variables in the frame of the center of mass defined at 
the required IPN accuracy. The individual positions y^ of the particles in the center- 
of-mass frame are given in terms of their relative position yi2 = yi — y2 and read at 
IPN order [25] 



yi 



ri2 



yi2 + Oic 



(2.6) 



together with 1^2 for the other particle. We introduced the total mass M = mi +m2, 
the relative mass ratios Xa = ruA/M, and the symmetric mass ratio u = X1X2 = 
mim2/M^, such that < ^ |, with z/ = | for an equal-mass binary, and z/ — )■ 
in the test-particle limit. We denote also the mass difference by 6M = mi — m2, and 
the relative mass difference by 6X = Xi — X2 = 6M/M, which can also be written 



in terms of the symmetric mass ratio as 5X 
velocity is V12 = dyi2/dt = Vi — V2, and f^2 = (vuVu) 



4z/. As previously the relative 
Note that the IPN correction 



in (2.6) vanishes for circular orbits. Similarly, from the time derivatives of (2.6) we get 



the IPN-accurate expression of the individual velocities as 



Vl 



Xo 



V 

2^ 



bX 



'12 



GM \ 
ri2 J 



V G5M 



V12 



2c2 



ni2^^i2)yi2 + 0{c 



(2.7) 



12 



The expressions (2.6)-(2.7) were derived in harmonic gauge, and we have still to 



check that they are also vahd in the new gauge specified by (2.3a)-(2.3&)- The center 



of-mass frame is defined by the vanishing of the center of mass position G. A generic 
gauge transformation will displace the position of the center of mass by the amount 

= — mi^i — m2^2, where ^a is the spatial gauge vector evaluated at the location of 
the particle A. The associated shift of the particle's individual positions with respect to 
the center of mass will then be given by ^^y^i = —5^G/M (the same for both particles). 
In the case of the gauge vector (2.3a)-(2.3&) we readily find ^1 = Gm2 1112/0^ and 
^2 = —6*^11112/0^, so that (5^G = and thus ^^y^i = 0. (See more details in the 
Appendix B of Ref. [26j.) 

We are now ready to write down the expansion of ta = |x — y^| when the CL 
ratio ri2/r tends to zero, where ri2 = |yi — y2| is the binary separation and r = |x| is 
the distance of the field point to the center of mass properly defined at the IPN order. 



Introducing the Legendre polynomials Pk and using (2.6), we thus have 

+00 



1 

ri 



fc=0 



1 + kXi 5X 



2c2 



X. 



r J 



Pfc(nni2) + 0(G,c- 



(2i 



together with 1^2. Here the argument of the Legendre polynomial is [nnu), the 
scalar product between n = x/r and ni2 = yi2/''"i2- Notice that the term proportional 
to G in Eq. (2.6) has consistently been neglected here. In addition, we see from 
Eqs. (2.5a)-(2.5c) that the scalar products (uaVa) are only needed with Newtonian 
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accuracy. Making use of (2.6)-( [278| , we get 



(riivi) = X2 



■ - ■ k 



[nvi2 



X2'y{nuvi2)] (^2^) Pkinnu) + 0{G,c-') . (2.9) 

fc=0 



We have a similar expansion for the term n\n-[/ri appearing in the spatial metric g, 



PN 
ij ■ 



Combining those results, we obtain the CL expansion of the 2PN metric (2.5a) 



(2.5c). In the companion paper p3j we shall use this expansion to compute the 
gravitational recoil generated during the ringdown phase. The gravitational recoil 
dominantly results from the coupling between the £ = 2 and i = 3 multipole moments 
[30] , where i is the azimuthal number of the Schwarzschild perturbation. Thus we need 
to expand the 2PN metric at least up to octupolar order, i.e. up to /c = 3. In fact it 
will turn out that it is necessary to push the CL expansion up to = 5 if we want to 
control all the terms which are of order the multipolar coefficients of the black 

hole perturbation in the Regge- Wheeler gauge. We shall discuss this point further in 
SecO 

At the zero-th order in the CL expansion we evidently recover the Schwarzschild 
metric of a black hole with mass M = mi + m2- Thanks to our gauge transformation 



(2.3a)-(2.3&) we find it to be directly given in usual Schwarzschild coordinates, namely 



2GM 



^7o'f™ = 0, (2.106) 



2GM 



^bchw ^ ^ -^n'n^ + 0{G') , (2.10c) 

where n = x/r is the unit vector in the direction of the field point x. As expected, 
we find that the Schwarzschild metric is exactly recovered in the limit z/ — t- 0, i.e. if 
one of the masses vanishes. Thus, we are now in a position to identify the 2PN metric 
expanded in the CL approximation with a perturbed Schwarzschild black hole with mass 
M, namely 

C = 9ir + Ku- (2.11) 
We find that the metric perturbation /i^i, expanded up to octupolar order in the CL 
approximation reads explicitly as 

^00 = ^^^[^{^^12) - 1] - v^^^{nni2Mnni2) - 3] + 3v—[2v^2 ' H12) ] 



G5M ri2 / \r 2 

-z/— 6{nvi2){ni2Vi2) + (^^12) 5fi2 - Q{nvi2) 

c^r r \ J 

+3z/(l - 3z/)^^ (t;?2[3(^^^i2)' - 1] - \{nvi2f\^{nni2f - 1] 

-("-12^^12)^ + '^{ni2Vx2){nvx2){nnx2) 
-^''(2-5<')^^-?.|5(™,,)=-31 
+9j,(1 _ 2i/)^^^2((™i2)(niji.i2)= + ^(nni2)(TO.i2)2[7(™i2)' - 3] 
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-Hi2)(^i2^i2)[5(nni2)' - 1]) + 0(G^c-^r^2) , (2.12a) 



-2^1 - 3z/)^^(nni2)[5(nni2)' - 3]t;i2 + 2z/^[t;?2 - Hi2)>i2 
-2z/(l - 3z/)^^— f2(nt;i2)(ni2fi2) + (n?2i2)[2fi2 - 3(nf i2)^])t;l2 



+1/(1 -2z/)—^ ^(2(ni2t;i2)[6(nfi2)(nr2i2) - (ni2fi2)] 



12 



+j/(2 - 7z/) — ^t;2^[3(nni2)' - \\v\^ 

-2v{\ - 7z/ + 13z/2)^^t;22[5(rmi2)' - 3]t;i2 
GM / 

+z/(l - 5z/ + 5z/^)^^ (^6(nni2)(ni2t;i2)^ + 5(r2ni2)Hi2)^[7(nr2i2)^ - 3] 



-6Hi2)(ni2t;i2)[5(nni2)' - 1] )t;l2 + 0{G\ c-\r\^ , (2.126) 
GMr^o/ . . . a .^ 3, 



/lij = 2z/-^^(^ni2n{2 - 6(nni2) n^a'^^) + -[5(n?2i2)^ - l]?2Vj 

GSMrfof, \ 7 i \2 T (« ii 

-6z/^ 5- (nni2) ^12^12 - [5(nni2) - ljn\2^-^^ 



5 

+ ;7('^'"-i2)[7(nni2)^ - 3]n*n-'' 
o 



GM 



c^r r 



^2[?,{nn,2f-l]v\2v{2-l 

+ {f^i2Vi2W{nni2){nvi2) - (ni2f 12)]'^*^ 



12 



- 3z.)^-f (^2[3(nni2)' - I]t;i2^l2 - i^{nv^2)Mnn^^^ - l^^ 



G6M r3 
5 



+1/(1 - 2z/)^^^(^-2(nni2)[5(nni2)' - 3]v\X2 



+ ^{nnu){nvuy[7{nnuy -3]d'' 

+3{nni2){ni2Vi2)^6'' - 3{nvu)ini2Vi2)[5{nnuf - IW' 

+9i;^2('^^i2)^i2^i2 ~ 9f ^2 [5(^2^12)^ — l]n[2'^'''' 

+ y^i'2(^m2)[7(nni2)'-3]nV) +0(G^c-^r^2). (2.12c) 

Even though we have performed the CL expansion including all terms up to order r^2 5 
only give in (2.12a)-(2.12c) the result up to octupolar order because of the proliferation 
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of terms at higher orders. But we further stress that those terms proportional to and 
r^2 in Eqs. (2.12a)-(2.12c) are fully under control (in an algebraic computer program), 
and were needed and used to get the final results given by Eqs. (3.6a)-(3.6/) and (3.7a)- 
( |3T/| ) below. 

If we come back for a moment to the three dimensionless scales epu, ^pn and ^cl 
defined in the Introduction, we can check that indeed the metric perturbation h^i, admits 
the general structure given by Eq. (1.4), in which the dimensionless coefficients /i^"'^'* 
are only functions of angles and mass ratios. 

Although the identification (2.11) we are making is mathematically cristal clear, 
we recall that its physical justification is not completely straightforward. We invoke 
the theory of the matching of asymptotic series, but use it in a formal way, since, as 
commented in the Introduction, the overlapping region between the domains of validity 
of the PN and CL expansions does not exist. Physically, we also rely on the fact that 
the merger as observed in numerical simulations lasts a very short time, which makes us 
feeling that the physics is essentially "conserved" when going from a PN description of 
the system to a perturbation of the final black hole. In addition, the PN approximation 
has proved to be very powerful in several past studies, with a domain of validity which 
often turned out to be larger than the one expected from elementary estimates (see 
pS] for a recent example in the extreme mass ratio regime). Here we are assuming 
a rather extreme extension of the domain of validity of the PN expansion — one for 
which > GM/c^, corresponding to the ultra-relativistic limit ^ 1. Nevertheless 
we shall find below and in the application [21] that the PN approximation performs 
well. 

In the next Section we shall describe our perturbation using the usual black 
hole perturbation formalism; for this we transform the coordinates from Cartesian 

= {ct,x,y,z} to spherical = {ct,r,6,(p}. This is appropriate for the 
spherically symmetric background, and the spherical coordinates are identified with 
the Schwarzschild(-Droste) coordinate system. Thus, 

^"^(^') = 9^5^^^-(^')- (2.13) 
We then write the scalar products (nni2) and (nw 12) in terms of them. Our conventions 
regarding orientations and various angles are explained in Fig. [Tj The unit vector 
n = x/r pointing from the center of mass to the field point, and the unit direction 
ni2 = 7x21^12 from body 2 to body 1 read 

n = (sin 6' cos v?, sin 6* sin yj, cos 6') , (2-14) 
ni2 = (cos/3,sin/3,0) . (2.15) 

Note our unconventional notation for the orbital phase angle /3. For a generic non- 
circular orbit, we have (see Fig. [T| 

(^^12) = sin 6* cos ((/? — /3) , (2.16a) 
{nvu) = sin 9 [fi2 cos (</}-/?)+ r^uj^ sin {(p - /?)] . (2.166) 
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Here the relative angular velocity is uJi2 = (3 with /3 being the orbital phase, and 
= (?T-i2^^i2) is the inspiral rate, where a dot stands for a derivative with respect 
to coordinate time t. Finally we find that all the components of the perturbation hai3 of 
the Schwarzschild metric 5'^^'^™ (both written in spherical coordinates x° = {ct, r, 6, (p}) 
are given as explicit functions of the spherical coordinates {r, 9, ip}, and depend on time 
t through the orbital parameters (3, Ui2, ri2 and ru- 



3. The 2PN metric in Regge-Wheeler-Zerilli formalism 

3.1. Multipole decomposition of a Schwarzschild perturbation 

We briefiy remind (see e.g. [3ll [32], |33l [3l] for more details) the usual decomposition 
into multipoles of a first-order perturbation of a Schwarzschild black hole of mass M. 
As usual, we write the perturbation as the sum of two kinds of perturbations, 

Kp = h^:l+h^:l (3.1) 

where the even-parity perturbation h^^^ essentially describes a perturbation along an 
(arbitrary) axis of the spherically symmetric Schwarzschild background, and where the 
odd-parity perturbation essentially describes a perturbation around that axisj^ Both 
perturbations are expanded with respect to a set of 10 tensorial spherical harmonics 
(cf. Appendix A). 

Following Regge and Wheeler's [35] conventions, the even-parity perturbation 
multipole decomposition reads |3l] 

4o^ = (i-^)E^o""^^-, (3.2a) 
hl^J =J2H'rYi,m, (3.26) 

£,m 



^oe = XI ^o^^eYi^m , (3.2c) 
/^g = E ^0 '"^^^^m , (3.26?) 



h':} = [i-—YY.H'ry^^m. (3.2e) 



i,m 



£,m 

h^:i = Y.h'rd^Y,,m, (3.2^) 

tt The even-parity perturbation is often called the "polar" perturbation, while the odd-parity one is 
called the "axial" perturbation. (From now on we pose G = c = 1.) 
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'hip — ' 



hie) 2 



l.m 

^ [ir^'" sin^ Q + G^''" (9J + sin ^ cos Q de)\ Y^^ 



(3.2z) 
(3.2j) 



where the summations over the integers i and m range from 2 to infinity, and from —I 
to I respectively. Note that the low multipoles £ = and £ = 1 correspond to the non- 
radiating pieces of the perturbation /i^^, and are not relevant to gravitational waves. 
For example, a monopolar perturbation [i = 0) would correspond to an infinitesimal 
shift of the black hole mass. (See e.g. [36t 132] for more details.) Similarly, the multipole 
decomposition of the odd-parity perturbation is 



h 



(o) 



'hp 



're 



fco sin 9 deYi^ 

i,m 



sin 9 

h'f:i = Y,k'r^^^9deY,^^, 

^2 ""^-^ {de - cot 9) d^Yi^m , 



i,m 



" sin 9 {de - cot 9) d^Y^m ■ 



(3.3a) 
(3.36) 

(3.3c) 
(3.3(i) 

(3.3e) 

(3.3/) 
(3.3^) 



Our convention for the scalar spherical harmonics Yg^rn is given in Eq. ( A.16 ) of Appendix 
In Eqs. ^ a[)-( |3.2j[ ) and ( |3.3a[ )-( |3.3g'| ) all the multipolar coefficients Hq"^, H^™', 
H2^, K^'"^, G^'™^ ^o'™5 ^i'*") ^o'*"' ^i"^ functions of {t,T} in Schwarzschild 

coordinates, and are defined in an arbitrary perturbative gauge. 

3.2. Computation of the multipole contributions 

Given the metric perturbation haB obtained from the CL approximation in the previous 



Section, we can compute all the coefficients ifg™, iff™. 



ki'^ . After some calculations 



consisting mostly of projections over the Zerilli-Mathews tensor spherical harmonics 



(using their orthonormality properties recalled in Appendix A), we obtain the even and 



odd multipolar coefficients in a particular gauge, which follows from the choice of gauge 



made in Eq. ( |2.3a[ )-( |2.36[ ). We shall however change to the Regge- Wheeler gauge, where 

'6 



the multipolar coefficients G^'™, /iq™", h'^{^ and ^2"^ vanish; this makes the expressions 



of the Regge- Wheeler and Zerilli master functions much simpler [see Eqs. (5.1a)-(5.16) 



below]. Note the difference between the choice of gauge (2.3a)-(2.36), which was made 
for the PN metric before its CL expansion, and a choice of gauge within black hole 
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perturbation theory, once the PN metric is in CL form. We transform the results to 
the Regge-Wheeler gauge by making the substitutions (see e.g. [Ml |37] for general 
expressions) : 



I'm 



H, 



HI 



£,m 




h: 



l,in 
2 



H: 



K 

£,m, 



£,m 
2 

£,m 







K 

'T£,m 

Kq 

£,m 



dth^ 



k\ 



£,Tn 



£,m 



(3.4a) 
(3.46) 
(3.4c) 
(3.4(i) 
(3.4e) 
(3.4/) 



We can now understand why it was necessary to expand the 2PN metric so as to include 



terms of order 



rj'2 in the initial perturbation (2.12a)-(2.12c). In the gauge 



transformation (3.4a)-(3.4/), the partial time derivatives of the multipolar coefficients 



in the initial gauge yield lower order powers of ri2 in the expression of the new multipolar 



coefficients in Regge- Wheeler's gauge. For example, from (3.4a) we observe that terms 
like G^'^^ ~ fi2/c^ ill the initial gauge produce terms like c~'^d^G^'^^ ~ 



m 



the multipolar coefficients H^'^^ in the Regge-Wheeler gauge. Such contributions of 
order ^^^d in the multipolar coefficients describing the perturbation in the initial 
gauge need to be consistently included in order to control all terms of order the 
multipolar coefficients describing the final perturbation in the Regge-Wheeler gauge. 

As expected, we find that the simplest of the Einstein field equations in vacuum in 
the Regge-Wheeler gauge is satisfied [up to terms (9(6"^, c~^, r^g)]; namely 



j^£,m _ j^£,m 



(3.5) 

Finally, we give below all the non-zero multipolar coefficients in the Regge-Wheeler 
gauge for all {£, m) up to £ = 3. All equations below are valid modulo remainder terms 
(9(G^, c~^, r^2)- For the even-parity perturbation. 



H 



2,0 



2,±2 



TT M 
V 

5 r 

677 M 
V 

5 r 




12 



12 



1 + (1 - 3z/) 



1 + (1 -3z/) 



9 

14 
9 



^^2 + 



17 

5 

14 
10, 



12^^12 



+ ^12 



14 



^"12 =F 777! '"12^^12^12 + 



+ (ri2Tiri2a;i2)'!>e^'''^ 



E 



3,±1 




371 5Mr 
V - 



12 



' 12 



19 



21 



1^12 T 3(1 



:^12'^12 



)] 



12^X2 



+ r 



2z/) iri2a;i2ri2 
1 



12^12 



3,±3 



± 




57r bMr 
V - 



12 f r 



12 



1 + 



12 =F g 1^2^12^12 + -r^^u 
19 



^^12^ (1 -2z/)iri2a;i2ri2 



-g(l + 7z/)rM2 



M 



1 o 



+ (ri2Tiri2a;i2)'U^'''' 



17 



(3.6a) 



(3.6&) 



(3.6c) 



(3.6d) 
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K 



■2,±2 



V — 
5 r 



K 



3,±1 



{^[l + (l-3-) 
-(rl2Tirl2a;l2)'}e^2^^ 



.2 ^ 10. . ^ 1 2 2 



[Sn 8Mri2(rl^\, 5/^ 19 \ .2 1/. o^• 



+ 



23 



, 2 2 



12"^12 



2 1 

''L T g i?^i2C^i2ri2 + gr^2a;?2 



K 



3,±3 



± 




^ 



5 / 19 \ o 
1 + g (^1 - y j T (1 - 2i/) 1 ri2a;i2ri2 



-^(1 + lv)rWi^ - (ri2 T iri2a;i2)'}e^3''^ , 



2,0 



^ y Mri2 . 
4W-t^ — 5— ri2 



9 

14' 



l + -(l-3z/K2 



(3.6y) 
(3.65) 

(3.6/i) 
(3.6z) 



' 67r Mri2 , x [ / \ 9 2 5 

— ^ (^2 =F 1 ri2UJi2) 1 + (1 - 3i/) — ri2 =F 



17 2 2 

+ |^^12^12 



14 



21 



1^12^12^12 



^1 ' = ±2y y [ri2 T - 1 ri2a;i2 + ^^12 ( g (^1 - y j ^2 



=F-(1 - 3z/)iri2Cc;i2 



(3.6i) 



(3.6^) 




1/ 



12 



(ri2 T iri2u;i2) 



6 



19 

y 



z/ r 



2 1 \ " 

=Fg(l - 2i/)iri2C^i2ri2 + g(l - llp)rl2Ujl2j ' 



12 

=F3i/3 



, g,- , (3.6/) 

where we made use of the relation — ^12 + '"12^12 when convenient. Similarly, the 
non-zero multipolar coefficients for the odd-parity perturbation are 



±2i 



yr 5Mrl 



12 



15 r 

13/ 68 

H 1 

28 V 13 



ri2'^12 



1 + 



9 



14 



5 

^12 T ~ 2i/)iri2<^i2ri2 



'3,0 



12 



ri2<^i2 



ON /5 37 73 2\ 



1^12 



/II 91 
^Vl8 ~ 18 



199 



2 / ,2 

12'^12 



r3,±2 



1 1 15% Mr^ 



14 



z/ 



12 



^12^12 



15 



V iri2a;i2ri2 + 



5 37 

6 ~ y 

3 



73 



(i-3i/) + \ ^-^y+^y~]r^2 



10 



6 

7 19 

2^ + y^") ^12^^12 



,2 , .2 



,=F2i^ 



V io r 



(3.7a) 

(3.76) 

(3.7c) 
(3.7(i) 
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ITC 



V 



3,±2 



^ — ^(1 - 3z/)(ri2 =F 1 ri2Wi2)ri2CUi2 e+ ^ 
4/ r"' 



(3.7e) 
(3.7/) 



One can check that any given multipolar coefficient F^'"^ in (3.6a)-(3.6/) and (3.7a)- 
(3.7/) satisfies the property F^"™- : 



where the overbar denotes the 



complex conjugation, consistently with the fact that the initial perturbation is real- 
valued. Furthermore, all the multipolar coefficients associated with the odd (or axial) 
perturbation vanish as expected in the zero angular momentum limit uji2 = 0, which 
corresponds to purely radial infall. We shall consider such head-on collisions in Sec. |5| 



4. Verification of the Einstein field equations 

As an important check of the previous results, we now verify that all the perturbative 
Einstein equations are satisfied. This requires the computation of the partial time 



derivatives of the multipolar coefficients (3.6a)-(3.6/) and (3.7a)-(3.7/). Recall that a 
generic multipolar coefficient F^'*" is function of the coordinate time t through the orbital 
phase f3(t), the orbital frequency (X'i2(t) = $(t), the distance ri2(t), and the inspiral rate 
ri2(t). Therefore, we have 



dtF 



dF 



Wl2- 



+ Wi2- 



dF 



+ ri2- 



dF 



+ ri2- 



dF 



£,m 



(4.1) 



dl3 dui2 dri2 dr^ 

The relative position, velocity, and acceleration of the two bodies can be expressed as 

yi2 = ri2ni2, (4.2a) 
vi2 = ri2ni2 + ri2u;i2Ai2, (4.26) 
ai2 = (ri2 - ri2 CU12) ni2 + (ri2 W12 + 2ri2 W12) A12 , (4.2c) 

where we have introduced the Frenet frame (ni2, A12) defined by A12 = L x ni2, with L 
being the unit vector orthogonal to the orbital plane, and in the same direction as the 
orbital angular momentum. Since we are working at linear order in G, the acceleration 
ai2 which is proportional to G can be neglected here, and we have 

ul,ru + OiG) , (4.3a) 
ri2 . .^.^^ ^^ 3^^ 



ri2 
^12 



-2ui2— + OiG). 
ri2 



Introducing these expressions of fi2 and UJ12 into (4.1), and neglecting terms 
0{G'^,c~^,rf2), we find for the even perturbation the non-zero partial time derivatives 



2,0 



dtH 



dtH^'^^ = 2 




ri2 



12 



(4.4a) 



^ ^ ^(^^12 =F l'^12Wl2) 

5 r-^ 



14 



^2 
12 



=F ^ir^uuru 



17 2 2 

+ ^^12^12 



(4.46) 
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3,±1 




7 r 
3 



2 \ J.2 



, 5/ 19 
(3ri2 =F 1 ri2a;i2) + 2 y ^ ^ ^12^12 



= ± 



^-(1 - 3i/)T;?2in2i^i2] +'yi2(^i2Tiri2i^i2)}e^'^, (4.4c) 

/stt 5m , X r ?^?2 r 5 / 19 \ 2 

V y 7^ (^12 T 1 ri2C^i2) 1 7^ [3 + - (^1 - -1/ j ri2 

^2(1 - 2z/) iri2a;i2ri2 + ^(1 - lluyl^u^ + {fi2 T iri2a;i2)'}e^''''(4.4cO 

(4.4e) 

/ N / 9 2 5 
1 + (1 - 3i/) (^Y^'^12 T ^iri2u;i2ri2 

(4.4i) 



= -4 W - i/ — — ri2 



1 + ^(1-3^^^2 1, 



Gtt Mri2 




(ri2 T i ri2u;i2) 



17 2 2 

+ |^^12^12 



Stt 5M (r 




j,2 



12 



(3fi2 T in2'^i2) + 2 ' ^ 



19 



2/ ^12^2 



9^^!^'° = 4 
9(^1'^^ = -2 

r3,±l 



^-(1 - 3i/)T;i2iri2u;i2] - ^;i2(ri2 =F in2^^i2)}e'^''' , 

/57r 5M .. . , f r^2 To 5 19 \ .0 
V y 7^ ^^12 T 1 ri2a;i2) ( 7^ [3 + - (^1 - -1/ j ri2 



(4.4p) 



=F2(1 - 2v) iri2C<;i2ri2 + 77(1 - 11^^)^126^ 



2 , .2 
12 



TT M 2 
5-^-12 

67r M 



l + ^(l-3i^K 



(ri2Tiri2a;i2)'}e^3^'^(4.4/i) 

(4.4i) 




(ri2 T i ri2C^i2)' [1 + - 3^)^?2] e^^'^ 




37r 5Mri2.. . . 
±A\l — v ^3 (ri2 Tin2'^i2) 



ri2 ± g i ri2u;i2 



+^?2 



5. _19 
6 V 5 
Stt 5Mri2 



1 



1 - ^i/)fi2 ± ^(1 - 5i/)iri2a;i2| e 



>/ 7 r'^ 




(ri2 =F 1^2^12)' 



5/ 19 

- 1 u |r 

6 V 5 



T^(l - 2i/)iri2a;i2ri2 + ^(1 - 5i/)rM2)] e^"''^ , 



and for the odd perturbation 



15 



9 13 \ 
IH 1 i/ 1 

14 V 3 / 



12 



-2 



TT Mr?, 



1/ 



12 



ri2'^12ri2 



/. ON /5 37 73 2\ 2 1 
(l-^-) +16 



3,±2 



2 /IStt Afr^ r 

^;^ri2a;i2(ri2 T i?^i2a;i2) (1 - 3z/) + ( ^ - ^z/ + 



14 r-' 

4 r r 2N. ■ /17 29 

^7^(1 - 5i/ + 5i/ )iri2a;i2ri2 + - 
15 V3U D 



5 37 

6 ~ "6 
65 2 
6 



73 



2 \ -2 

z/ )ri2 



Z/+ — z/ )r^2^^i2 



2 , ,2 



==F2i/3 



(4.4j) 
(4.4A;) 

(4.40 

(4.5 a) 
(4.56) 

(4.5c) 
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dtK 
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3,0 



-(1 - 3z/)f 12^^12^12 , 



3,±2 



'Svr Mri2 



1 - 3z/)(ri2 =F i'^i2Wi2)^n2Wi2 e'F^''' 



19 

(4.5e) 



42 r2 

We then check that the Einstein equations are satisfied for all (£, m) up to £ = 3 
for a generic non circular orbit, up to terms c~^, r^2)- We give them here in the 

linear case for completeness (see e.g. [361 EH EH] for general expressions). For the even 
perturbation, these seven equations read 



Stiff'" - driH''"' - K 







'r£,m 



dt 
dt 



2r2 ^1 



r 



2r2 







+ 1) 



+ 1 



H 



dt 



0, 



r 



Lm 



K 







(4.6a) 
(4.66) 

(4.6c) 



(4.6(i) 



(4.6e) 



(4.6/) 



together with Eq. (3.5). Note that these equations are not all independent. For example, 
if Eqs. (4.6a) and (4.66) are satisfied, then Eq. (4.6/) is also satisfied. For the odd 
perturbation, the three remaining equations are 







dt 
dt 



r 



^2 '^O 







dtK 



(£-!)(£ + 2) -,, 



0. 



(4.7a) 
(4.76) 

(4.7c) 



Assuming that all terms of a given equation in (4.6a)-(4.6/) and (4.7a)-(4.7c) are of the 
same order of magnitude, we can now understand by coming back to the initial metric 
decomposition (3.2a)-(3.2j) and (3.3a)-(3.3g') that it is necessary to expand 5^™, g™ 



and g™ at the same PN order in Eqs. (2.5a)-(2.5c), i.e. up to 0{c~^) included. The 
previous verification of the field equations provides a good check of the algebra yielding 
the perturbation coefficients (3.6a)-(3.6/) and (3.7a)-(3.7/). 
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5. Numerical evolution of the perturbation 
5.1. Regge- Wheeler and Zerilli master functions 

From the multipolar coefficients if^'™, iff™, K^'"^, /cq™ and kf"^ one can construct for 
any {i, m) two gauge-invariant scalar fields, namely the Regge- Wheeler [55j function 



^'^^ and the Zerilli [36] function ^I/^^, which contain all the information about the 
perturbation of the Schwarzschild metric. Gauge-invariant expressions of '^f^'' in terms 
of the multipolar coefficients in a general gauge are given e.g. in [321 |33l [39]. In the 
Regge- Wheeler gauge [35], the coefficients G^'"^, h^"^, and ^'"^ vanish, so that these 
expressions get simplified and read 

r - 2M 



(5.1a) 



where we introduced the widely used notation = |(£ — 1)(£ + 2) [ID]. Note that 
the multipolar coefficent H^'"^ does not enter the expression of '^fl^. Because we are 
considering linear perturbations, the master functions "^f^^ are defined up to a scale 
factor. We use the same convention as in ^3], emphasizing the link between "^f^^ and 
the polarization states and of the gravitational waves at future null infinity; with 
our convention the two independent -|- and x polarization states are given by 



^+ - i = ^ E y 1^ + i + ^(^") ' (5-2) 

where -2^^.m denotes the spin-weighted spherical harmonics of weight —2. The 
asymptotic waveform is also related to the more fundamental Weyl scalar \l/4, which 
admits a closed-form expression in terms of the master functions "^f^^ (see the 
Appendix). 

The two master functions satisfy a wave equation with specific potentials V^^'"'', 

{d^-dl + vt''^)¥;-^ = o, (5.3) 

where the so-called tortoise coordinate is related to the Schwarzschild radial 
coordinate r by 

r, =r + 2Mln(^-l). (5.4) 
Notice that these wave equations are only valid in vacuum; otherwise one has to include 



a source term in the right-hand-side of (5.3), see e.g. [32l |33]. The Zerilli and Regge- 
Wheeler potentials read respectively 

(5.5) 
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with 



(e) 



Xe{Xe + 2y + 3M(r - M) 



(5.6a) 



(o) 



1. (5.66) 

One can easily prove that | < < 2 for all £ ^ 2 and for all r such that 

2M < r < +00, showing that the potentials V^'^^ and V^"^ are very similar |10j. 



5.2. Numerical evolution 



The wave equations (|5.3|) are evolved numerically using the initial conditions at time 



t = and for any tortoise radius r*, namely ^^'^^lO' ''^*) ^^"^ ^t'^fmi^^ derived from 
the CL expansion of the 2PN metric for compact (i.e. point-mass) binaries. These 



initial conditions are calculated by plugging Eqs. (3.6a)-(3.6/) and (3.7a)-(3.7J) and 



their partial time derivatives (4.4a)-(4.4l) and (4.5a)-(4.5e) into (5.1a)-(5.16) and their 
partial time derivatives. 



We use Dirichlet boundary conditions, setting "^f^^ {t, rl 



(c,o). 



at some radii r™''^ and rl 



We choose the radii r™™ and 



If J. „max\ _ n 
'£,m I''' ' * ) — ^ 

in such a way that 



these boundary conditions are causaly disconnected from the computational domain, 
i.e. the spurious radiation generated on the boundaries {t, 7-™°'™ax j ^Qgg ^q^^ have time 



to propagate up to the extraction radius 



ext 



1 /^max 

2 V * 



for ^ t ^ t™'^'', where 



Extending the computational domain to [r™ 



require using the Sommerfeld boundary conditions [H] {dt + dr^)'^f^\t,r 



would 



{dt — drj'^f^\t,r'^™) = 0, which are approximate boundary conditions, or even better 



some exact boundary conditions |l2l |13] . The results below are based on computations 
where we have chosen r™'"^ = — 60M and r™*™ = 660M, such that rf^^ = 300M and 
^ma^ = 360M. 

A simple explicit second-order finite difference scheme has been used to evolve 
the wave equations (5.3). We always choose the spatial grid resolution 6r^ and the 



time increment 6t such that the so-called Courant-Friedrichs-Lewy condition 6t < 6r^ is 
verified; therefore the code is stable. The results in Figs |3||4| [5] and |6] below are based on 
computations where we used a spatial grid resolution (5r* = 0.2M, and a time increment 
5t = O.IM. 

In order to check the second-order convergence of the code, we computed (for 
example) the real part of the {£, m) = (2, 2) mode of the Zerilli master function, 
i/j = ^[^'^\], for different spatial grid resolutions Sr^, = 0.2M/h, where h = 1,2,4, 
with a constant time increment St = 0.025M, in the case of an unequal mass binary 
on circular orbit with u = 0.185, ri2 = 1.6M and /3 = 0. The good overlapping of 
the differences ip\h=i — 4'\h=2 and 4:{ip\h=2 — 4'\h=i) as shown in Fig. [2] demonstrates the 
second-order accuracy of the code. 

The numerical code was tested in several ways, and against previous published work 
as well as on some material presented in [21]: 



Close-limit approximation with post- Newtonian initial conditions 



22 



0.00015 



0.0001 



5e-05 



-0.0002 



-0.00025 



-0.0003 



280 



¥l/,=i -¥l/,=2 

4(¥lA=2-¥lft=4) 




290 



300 



310 320 

? (units of M) 



330 



340 



350 



Figure 2. The good overlaping of the differences ip\h=i — i'\h=2 and A{%}j\h=2 — V'U=4) 
demonstrates the second-order accuracy of the code. 



(i) 



111 
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Using the Misner initial data [44] as given by Price & Pullin [T7j, we reproduced 
in the case of head-on colhsion the waveform of their Fig. 2 and the associated 
radiated energy; cf. Eq. (16) in [T7]. 

Using the initial data provided by Sopuerta et al. [20] (a conformally flat 3-metric 
with a Bowen-York extrinsic curvature and a Brill-Lindquist conformal factor), we 
reproduced their waveforms in Fig. 7, and the fluxes of energy, angular momentum 
and linear momentum of their Figs. 4, 5 and 8 respectivelyjft] 
Checking that the total energy, angular momentum and linear momentum radiated 
do not depend on the physically irrelevant initial phase (3; and that the components 
of the integrated linear momentum flux, or gravitational recoil, transform according 
to the usual law for vectors under a shift of the initial phase /3. 

Checking that the quasi-normal mode frequencies of the waveforms are in good 
agreement with theoretical values [311 HSj US] . 



ft The updated plots of Fig. 7 and Fig. 8 of [20 (taking into account the corrections from their first 
Erratum) are available on the e-Print server arXiv.org. Fig. 7 actually shows the Zerilli-Moncrief and 
Cunningham-Price-Moncrief master functions, and not their time derivatives as stated. 
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5.3. Ringdown waveforms for head-on collisions 



We first consider head-on collisions for wliicli tlie perturbation is purely polar, i.e. 
^£°m ~ 0- thus show in Figs. [3]and|4]the real part of the Zerilli master function for 
even-parity (or polar) perturbations, Initial data for this case are obtained 

by setting Uu = (no orbital angular momentum) and ri2 = (time-symmetric 



initial conditions) in the expressions (3.6a)-(3.6/) and (4.4a)-(4.4/) of the even-parity 



multipolar coefficients and their partial time derivatives. 

In Fig. [3]we consider an equal mass binary (z/ = |), and compare our 2PN-accurate 
results to those of Abrahams & Price [H], who studied similar head-on collisions 
using Brill-Lindquist (BL) initial data |17]. We shall restrict the comparison to the 
{£,m) = (2,0) mode ^2^0 for simplicity. This comparison requires a detailed discussion 
of the relation between the two notions of distance between the two bodies used in both 
initial data sets. Using our conventions for the perturbation and various angles, we find 
that the only non-vanishing multipolar coefficients in the BL geometry are 

1 f^ML"^ 1 



2,0 

BL 



K 



2,0 

BL 



5 i?3 I 



M_ ' 
2R 



(5.7) 



where R = \ {y/r + y/r — 2M)^ is the isotropic radial coordinate, and L is the distance 
between the two black holes of the BL solution. Because the multipolar coefficients G 



Hq^ and /i^'™ vanish, the gauge in which the perturbation (5.7) is written coincides with 
the Regge- Wheeler gauge (hence our use of the symbol ~ on the multipolar coefficients). 
Setting u:i2 = and ru = in (3.6a) and (3.6e), we get for our PN initial conditions 

~2 



2,0 

PN 



K 



2,0 

PN 



TcMrj^ 



where we used r 



5 r3 2 V 5 i?3 ^1 + 

i? (1 + H)^. Observe first that if we set ri2 



2r) 



(5.8) 



L, then the two 



perturbations (5.7) and (5.8) coincide as they should in the weak-field domain R M. 



But in the strong field domain R > M/2, the comparison of the initial distances L and 
ri2 is difficult. It then becomes interesting to check if these two measures of the distance 
between the holes can be related in such a way that the two waveforms compare well. 
Most of the perturbation '^'"20 that propagates to future null infinity is generated 

: 2 potential V2°^ for polar perturbations, which is located 



around the maximum of the I 
around r ~ 3.1M, or in terms of the isotropic coordinate R 



2M. If we wish to identify 



the perturbations (5.7) and (5.8), it is then natural to impose the definition 



ri2 
L 



2r) 



5/2 



1.75. 



(5.9) 



R~2M 



We show in Fig. [3] the (£, m) = (2, 0) mode of the Zerilli master function using both 
BL initial data with L = 2M, and our PN initial conditions with ri2 = 1.75L. We 
observe that the waveforms compare very well, which means that our post-Newtonian 
initial conditions are essentially equivalent to the Brill-Lindquist initial data in the case 
of head-on collisions. We find that the waveform computed with PN initial conditions 




280 290 300 310 320 330 340 350 

? (units of M) 

Figure 3. The real part of the mode (2, 0) of the Zerilh master function for even-parity 
(or polar) pertmbations, in the case of an equal mass {i^ = j) head-on collision, using 
Brill-Lindquist initial data as given in Ref. [18] with an initial distance L = 2M (red), 
and our PN initial conditions with ri2 = 1.75L (green). 



is slightly delayed with respect to the BL one. So in order to achieve this agreement we 
also had to translate in time the PN curve by an amount At ~ 4M. We checked that 
this good agreement does not depend on the value of the initial distance L. 

We now focus on the waveforms obtained from our CL approximation with 2PN- 
accurate initial data. In Fig. |4] we consider the head-on collision (with time-symmetric 
initial conditions) of an unequal mass binary with mass ratio v = 0.2, and show the 
real part of the non-vanishing modes of the Zerilli master function for even-parity 
perturbations, 3f?[^^^m]; up to £ = 3. We set the initial distance to a value twice larger 
than in the orbital case shown later, i.e. ri2 = 4.4M, so that the Newtonian energy of 
the binary E = —Gmim2/ri2 is identical to that of the circular orbit configuration with 
initial distance ri2 = 2.2M. This configuration being axisymmetric, for a given i all 
{i, m) modes can be related to the (£, 0) mode. For example we know from Ref. [H] that 

= and = -^1^^. We checked that the ratios of the amplitudes of 

the modes shown in Fig. |4] are in very good agreement with these theoretical values. 

5.4- Ring down waveforms for circular orbits 



In the case of circular orbits, both even and odd perturbations contribute. In this case 
we simply have to set ri2 = in Eqs. (3.6a)-(3.6/), (3.7a)-(3.7/), (4.4a)-(4.4l) and 
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Figure 4. The real part of the non-vanishing modes {£, m) of the ZeriUi master 
function for even-parity (or polar) perturbations, up to £ = 3, in the case of a head-on 
collision with v = 0.2 and ri2 = 'iAM. 



(4.5a)-(4.5e). Our initial data will therefore depend only on the (physically irrelevant) 
initial orbital phase /3, the initial distance ri2 and the initial orbital frequency uji2- We 
know from the IPN-accurate equations of motion that, for a circular orbit, the orbital 
frequency is related to the binary's separation by the Kepler-like law [6] 



M 



12 



l + (z/-3) 



M 
ri2 



+ 0{c 



-4\ 



(5.10) 



so we have only one free parameter ri2. This initial orbital distance ri2 will be an 
important parameter since it will be used in applications like [21] as a ^''matching radius" 
to connect the computation of the ringdown phase to the previous inspiral and/or plunge 
phases. 

An important point is worth emphasizing at this stage. Recall that in our 
previous calculation of the initial data for the Regge- Wheeler and Zerilli equations we 
have systematically and consistently neglected the non-linear terms 0{G'^). Thus the 
perturbation coefficients (3.6a)-(3.6/) and (3.7a)-(3.7y) we consider, and which are valid 
for general orbits, are linear. Now, by introducing the expression (5.10) of the orbital 
frequency Ui2 (where M = GM) into the results (|3^-(|33|, ^J^ -^J]^, ( |CT| )-(|44]|) 
and (4.5a)-(4.5e) for the multipoles and their time derivatives, we do obtain terms which 
are of order 0{G^) or more in the case of circular orbits. Those terms have to be kept as 
they are, because the result (5.10) comes from an independent calculation at the level 
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Figure 5. The real part of the non- vanishing modes {£,m) of the ZeriUi master 
function for even-parity (or polar) perturbations, up to £ = 3, in the case of a circular 
orbit with j/ = 0.2, ri2 = 2.2M and 13 = 0. 



of the equations of motion. That is, once we have proved (in Sec. |4]) that the Einstein 
field equations are satisfied for generic non-circular orbits, we are allowed to reduce the 
solution to the particular case of a circular orbit by inserting (5.10); our point is that 
this adds new powers of G which constitute an integral part of our solution of the field 
equations. This being said, a more involved calculation making use of the theory of 
second-order perturbations of a Schwarzschild black hole would introduce other terms 
of the same order in the final solution. But we do not have access to these terms in 
this work, which is based on first-order perturbations. 

We show in Fig. |5] the real part of the non-vanishing modes of the Zerilli master 
function for even-parity perturbations, S^i^^''^]' up to £ = 3, in the case of the unequal 
mass binary with mass ratio u = 0.2 (same as before), on a circular orbit with initial 
separation ri2 = 2.2M and initial phase 13 = 0. Notice the peculiar behaviour of the 
i = 2 modes which do not vanish asymptotically. This is because some terms in the CL- 
expanded 2PN metric (2.12a)-(2.12c) behave as 1/r in the far zone, therefore yielding 
non- vanishing asymptotic values for the quadrupolar modes of the Zerilli master function 
[see e.g. the third term in hoo in Eq. (2.12a) . 

The Newtonian energy of the binary on a circular orbit in Fig. |5] is chosen to be the 
same as for the head-on collision shown in Fig. |4| Notice the difference in amplitude, 
which is typically a factor 2, between the orbital case and the head-on case; the even 
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Figure 6. The real part of the non-vanishing modes {£, m) of the Regge- Wheeler 

master function for odd-parity (or axial) perturbations, up to = 3, in the case of a 
circular orbit with v = 0.2, ri2 = 2.2M and /3 = 0. 

(polar) perturbation is stronger in the head-on colhsion case than in the orbital case 
for the same total energy. This is consistent with the fact that in the case of a circular 
orbit, we have also in addition to the even-parity perturbations, some odd-parity or 
axial perturbations. 

The odd/axial modes for the circular orbits are described by the Regge- Wheeler 
master function ^E'^"^. Notice that they can be computed with our PN initial conditions; 
Brill-Lindquist (BL) and Misner-Lindquist (ML) initial conditions are time-symmetric 
and can be applied only to the polar modes of head-on collisions. (See [201 [21] for initial 
conditions generalizing BL and ML, and which also permit to compute the axial modes 
of circular orbits in the CL approximation.) It is therefore particularly interesting to 
compute the real part of the Regge-Wheeler function, 3?[^I/^°^], as we do in Fig. |6] for the 
same unequal mass binary, and for the same initial conditions as for the polar modes 
shown in Fig. [Sj Comparing Figs. [5] and |6] we find that the amplitude of the even modes 
is approximately one order of magnitude larger than the amplitude of the odd modes. 
The even modes, which were the only ones present for head-on collisions without angular 
momentum, still dominate when we turn on the orbital angular momentum and consider 
the circular orbits. 
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6. Conclusion 

We have proposed an implementation of the close-limit (CL) approximation for binary 
black holes starting from post-Newtonian (PN) initial conditions developed at 2PN 
order. The 2PN metric generated by two point particles, when restricted to be linear 
in G in order to be consistent in fine with the linear black hole perturbation, was 
formally developed in CL form and identified with the metric of a linearly perturbed 
Schwarzschild black hole. We proved that the resulting multipolar coefficients describing 
the even-parity and odd-parity components of that perturbation consistently satisfy the 
perturbative Einstein field equations for general non-circular binary orbits. 

The post-Newtonian initial data was then specialized to the cases of head-on 
collisions and of circular orbits. We evolved numerically the Regge- Wheeler and Zerilli 
wave equations starting from those initial data, obtained the waveform generated during 
the ringdown phase, and compared these two cases. In a separate work |2l] , we apply this 
formalism to the computation of the gravitational recoil produced during the ringdown 
phase of circular-orbit compact binaries, and match it to a previous PN calculation of 
the recoil accumulated in the preceding inspiral and merger phases |49j . 

They are several ways in which this work could be extended. First, one may expand 
to a higher order in the CL approximation, to get higher multipoles of the perturbation, 
or include in the initial 2PN metric terms describing the effect of spins of the initial 
black holes. We could also include terms or higher in the initial post-Newtonian 

metric, but it would be necessary to use the theory of at least second-order perturbations 
of a Schwarzschild black hole |28j . 

An important limitation of the present calculation is the impossibility to treat 
perturbations of a highly spinning black hole. In this work the initial orbital angular 
momentum of the binary had to be included in the perturbation of the final non-rotating 
black hole. A solution would be to employ the Teukolsky formalism [50] to perform 
similar calculations using a Kerr background instead of a Schwarzschild one. This 
would allow for a better description of the ringdown phase of the final black hole which 
is known from numerical calculations to be a rapidly spinning Kerr black hole [71 El |9] . 
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Appendix A. Basic material for black hole perturbations 

Although the material contained in this Appendix is well-known, we give in self- 
contained form the usual decomposition of an arbitrary linear perturbation ha^ of a 



Close-limit approximation with post- Newtonian initial conditions 



29 



Schwarzschild black hole onto the Zerilli-Mathews tensor spherical harmonics (correcting 
misprints frequently found in the litterature) , and recall the relation between the Weyl 
scalar and the Regge- Wheeler and Zerilli master functions, from which we derive the 
asymptotic waveform in the form of its two polarization states + and x . 



Appendix A.l. The basis of tensorial spherical harmonics 

In order to perform practical calculations, it is particularly convenient, instead of using 



the decomposition (3.2a)-(3.2j) and (3.3a)-(3.3g') introduced in Sec. |3| to introduce the 
Zerilli-Mathews basis of tensor spherical harmonics {e^'^'"^)A=i,-,io, and to write the 
perturbation as 

10 



(A.i; 



A=l e,m 

where = {t, r, 9, ip} are the usual Schwarzschild spherical coordinates. This is because 
the tensor spherical harmonics (e^^^'™)^=i_... 10 form an orthonormal basis [511 ESI 
in the sense that 



I C^C' e^^-'" ef/'"'' dn = 5aa' ?>ram' , (A.2) 

where is the Euclidean metric written in spherical coordinates, bij is the usual 
Kronecker symbol, the overbar denotes complex conjugation, and the integration is 
performed on the sphere of unit radius. This basis being orthonormal, each component 
hA,i,m of a given perturbation hap can be calculated by the projection 



hA,e,r. 



^A,e,m 



h), 



from which we then deduce the multipolar coefficients Hq"^, iff™. 



(A.3) 

^2 defined by 



(A.17)-(A.26) below. This method was systematically used to get the results (3.6a)- 
(3.6/) and (3.7a)-(3.7/) from the information on the metric perturbation contained in 
Eqs. ( |2l2^ -( |2l2^ . 

The Zerilli-Mathews basis of tensorial harmonics explicitly reads 
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(A.7) 



(A. 



(A.9) 



(A.IO) 



(A.ll) 



(A. 12) 



(A.13) 



-SV2 -sVi J 

where we introduced the convenient shortcut s = sin 6*, together with the angular 
operators 



Vi = 2 {de - cot 9) d, 
Vo = dl - cot 6 de 



'2 o2 



d 



and where our convention regarding the spherical harmonics is 



Y, 



l2i+l {i-m)\ 
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(A.14) 
(A.15) 



(A.16) 



47r {£ + m)\ 

with P£m{x) = (— 1)™(1 — a;^)™/^(^)™'P^(x) being the associated Legendre functions. 
Finally, the coefficients hA,e,m of the arbitrary perturbation ha^ are related to 



the multipolar coefficients H^"^, H^'"^, 



h 



£,m 



of the decomposition (3.2a)-(3.2j) and 
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(3.3a)-(3.3g) through 

hi,i,m = — l/fo'''", (A. 17) 



2M\ , 



r 



h2,i,m =-iV2Hl'"', (A. 18) 

h3,e,m =(l-—) H'r. (A. 19) 



r 



/^4,,„. =-i^/2|±j|/.^'-, (A.20) 
h,,,,m =1/2^^^ /i'r, (A.21) 
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Appendix A. 2. Link between the asymptotic waveform and the master functions 

Here we recall the general expression (known as a Chandrasekhar transformation |10] ) 
of the Weyl scalar \l/4 in terms of the Regge- Wheeler and Zerilli master functions 
and we compute the combination h+ — ihx, where /i+^x denote the two asymptotic wave 
polarizations. 

In the Schwarzschild spherical coordinate system {t,r,9,ip}, we let (e",eg,e^) be 
the associated orthonormal basis of t = const hyper surf aces, and ef be the timelike unit 
vector orthogonal to them. Then, we consider the following complex null tetrad: fc" 
and Z° are two real null vectors, while m° and m" are complex conjugated null vectors 
defined by 

r = (er + e'^) , (A.27) 
r =-^(e^-e?), (A.28) 
m" = -^(e^-ie3, (A.29) 



32 
(A.30) 
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We therefore have rrf'ffia = 1 = —l^'ka, all the other scalar products vanishing. Various 
conventions for the definition of the Weyl scalar \E'4 can be found in the litterature. We 
adopt 

^4 = C^^^5 r r m^ (A.31) 

where Cap-ys is the Weyl tensor, which coincides with the Riemann tensor in vacuum. 
It can be shown that, for a generic perturbation of a Schwarzschild black hole. 
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(A.32) 



where r is the usual Schwarzschild radial coordinate, and the master functions 



(e,o) 

are those defined in Eqs. (5.1a)-(5.16). The relation (A.32) is exact for first-order 



perturbations of the Schwarzschild geometry [40]. The s = —2 spin- weighted spherical 
harmonics -2^,m. are defined for any integer s by [531 El] 
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where the Wigner functions gde^m read 



with fcmin 
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(A.34) 



max(0, m — s) and k„ 



= min(£ + m, £ — s). In the case s = we recover 
the scalar spherical harmonics (A. 16). The differential operators acting on "^^^^^ in the 



relation (A.32) read explicitly as 



(e,o) 



where r^, is the tortoise coordinate (|5.4|), and the potentials wf''"^ are given by 
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(A.35) 

(A.36) 
(A.37) 

Now, in the limit r — )■ +oo, the Regge- Wheeler and Zerilli functions "^f^^ are 
functions of th e retarded time t — r^, only. Indeed, they are solutions of the wave 
equations (5.3) with decaying potentials v]*^'"^ oc l/r^. Because the potentials W^'^'°^ 
also decay as r^^, we have the asymptotic expressions 



Recall that = ^ (£-!)(£ + 2). 
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(A.38) 
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Furthermore, working in the transverse and traceless gauge, and performing some 



projections onto the hnearized Weyl tensor around flat spacetime in the definition (A. 31 ) 
leads to the well-known expression 

^4 = 9,2(/i+-i/ix) + 0(r-'), (A.39) 
where the two polarization states /i+ and are defined by 

h+ = l {^e4 - e;ey hi^ = ^ {hee - V) > (A.40) 
hx = \ (44 + e^s) hij = hg^ . (A.41) 
Finally, combining the results ( A.32[ ) and (A.39) with the asymptotic expansion ( |A.38 ), 



we recover the well-known formula 
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